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ABSTRACT 


New  variational  bounds  are  derived  on  the  generalized  polarizabilities 
of  a quantum-mechanical  system,  for  arbitrary  complex  frequencies  c = v+iw 
and  two  different  perturbations  u and  v.  No  power  of  the  Hamiltonian  h 
higher  than  h^  is  involved  in  the  bounding  functionals.  For  a certain 
range  of  v-values,  upper  and  lower  bounding  functionals  are  obtained  which 
contain  merely  a single  trial  vector  but  also  introduce  an  inverse  operator  like 
h*^.  This  impractical  feature  can  be  avoided  with  a subsidiary  variational 
principle,  leading  to  bivariational  upper  and  lower  bounds.  Explicit 
bivariational  bounds  are  also  derived  which  are  valid  for  all  values  of  c . 

Both  theoretical  and  practical  aspects  of  the  bounds  are  discussed. 
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NEW  VARIATIONAL  BOUNDS  ON  GENERALIZED  POLARIZABILITIES^ 

Peter  D.  Robinson* 

1 . INTRODUCTION 

Let  a quantum-mechanical  system  be  described  by  a self-adjoint 
Hamiltonian  operator  h in  a complex  Hilbert  space  H , and  suppose  that  h 
possesses  a complete  set  of  orthonormal  eigenvectors  {6|^}  with 
corresponding  energy  eigenvalues  {Ej^}  . If  the  system  is  in  a state  0^, 
its  dynamic  polarizability  a(c)  at  complex  frequency  c = v+ioj 
associated  with  a perturbation  u can  be  defined  as 


a(0 


2 Real  I (E.-E  +0 
Mn  ^ " 


<uen.ek><ek.iien>  . 


(1.1) 


the  summation  being  over  all  states  different  from  6 . The  notation 

n 

<,>  denotes  the  complex  inner  product,  so  that  for  all  il>  and  4*  in  H 
and  complex  numbers  s we  have 


= <4',4»,  <S*,4'>  ■ S<4>,4'>,  «I>,S4'>  = S<3>,4'>,  (1.2) 

a bar  denoting  complex  conjugate.  Previous  authors  [see  for  example  references 
and  ^-27^]  have  presented  bounding  variational  functionals  on  a(v) 
or  a(iu))  , the  dynamic  polarizabilities  when  i,  is  wholly  real  or  wholly 
imaginary.  Often  there  has  been  a restriction  to  real  u and  real  6^^.  With 
shortcomings  of  many  of  these  bounding  functionals  have  been  the  high  powers 
of  h involved,  and  a multiplicative  factor  of  w'*  which  is  unfortunate  for 
small  u>. 

Expressions  similar  to  that  in  (1.1),  but  with  inner  products 

<v0n.9k^'^®l(*‘^®n^  involving  different  perturbations  u and  v,  define  quantities 
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arising  for  example  in  theories  of  optical  rotatary  power  [9],  and  nuclear- 
magnetic  shielding  and  chemical  shifts  [10, IJ].  Such  more  general 
expressions  can  also  arise  in  double  perturbation  theory  [1_2].  With 
generalizations  such  as  this  in  mind,  as  well  as  the  desirability  of 
admitting  arbitrary  complex  frequencies,  perturbations  and  unperturbed 
states,  we  show  how  to  derive  upper  and  lower  bounding  variational 
functionals  on  the  quantity 

Z(f,g;c)  = 2 Real  I (E.-E  <g,0^><e.  ,f>  . (1.3) 

k^n  K K 

No  separate  consideration  is  necessary  for  quantities  defined  as  imaginary 
parts  of  summations  like  that  in  (1.3),  for  g can  merely  be  replaced  by 
-ig  if  necessary.  Without  significant  loss  of  generality,  the  arbitrary 
complex  vectors  f and  g are  taken  as  members  of  the  reduced  Hilbert 
space  WpC  W , containing  all  vectors  in  H whichare  orthogonal  to 
6„  . I.e, 

f,g  € <*,9p>  = 0)  • (1.4) 

Apart  from  their  intrinsic  theoretical  interest,  bounding  variational 
functionals  can  in  principle  lead  (with  suitably  artificial  choice  of  trial 
vector)  to  bounds  on  unknown  quantities  in  terms  of  certain  known  quantities, 
like  sum  rules  or  moments  [4,6].  However,  if  they  are  to  be  a viable 
practical  tool,  bounding  functionals  must  not  present  exceptionally  severe 
problems  of  evaluation  when  reasonable  trial  vectors  are  employed.  It  is 
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for  this  reason  that  we  do  not  much  concern  ourselves  with  functionals 
which  involve  powers  of  the  operator  h higher  than  the  second. 

2.  A BIVARIATIONAL  APPROXIMATION  TO  Z 

Variational  approaches  to  the  task  of  bounding  Z stem  from  its 
alternative  but  equivalent  specification  in  terms  of  the  solution-vector 
(t>  of  the  equation  in 

{h-E^n)0>  = f , <^.f  c . (2.1) 


This  is  simply 

Z(f,g;^)  = 2 Real  <g,((»  = <q,(j)>  + <(|),g>  . (2.2) 

Setting  c = v+ iw  (with  v and  oi  real),  and  defining  for  convenience 


H » h-E^+v  , (2.3) 

equation  (2.1)  is 


A«(>  = f . ♦.f  c (2.4) 

with 

A = H + iu;  . (2.5) 
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t 


This  decomposition  of  the  linear  operator  A as  the  sum  of  a self-adjoint 
part  H and  a skew-self-adjoint  part  iw  is  important  for  the  establishment 
of  the  bounds  in  §3. 


Along  with  equation  (2.4)  we  consider  the  auxiliary  equation 


where 


A*<|^  = g , il'.g  G 

A*  = H - id) 


Is  the  Hilbert-space  adjoint  of  A.  We  note  that 
«|»,g>  = «|),A*i|/>  = <A4»,\|^>  = <f,i|;>  , 

SO  that  we  can  express  Z in  the  form 

Z(f,g;c)  = <g,'l»  + <f.'l^> 

in  terms  of  the  solution-vectors  4»  and  iii  of  equations  (2.4)  and  (2 
Associated  with  this  pair  of  equations  is  the  bivariational 
approximation  to  <g,^>  given  by 

R(1',*)  * -<'f,A*>  ♦ <’f,f>  <g,^>  , y,*  E , 

with  the  trial  vector  f playing  the  role  of  a kind  of  Lagrange 
multiplier.  In  terms  of  the  difference  vectors 


(2.6) 

(2.7) 

(2.8) 

(2.9) 
.6). 

(2.10) 
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6it>  = e 6(p  = <>-(}.  c , (2.11) 

we  have  the  relation 

R(4',<I')  = <g,((»  - <6tl^,A64»  . (2.12) 

The  complex  conjugate  of  R('f,'l))  is 

fi('i',<l>)  = -<4>,A*'}'>  + <4',g>  + <f,4'>  (2.13) 

= <f,V*>  - <6(p,A*S\p>  . (2.14) 

which  is  a bivariational  approximation  to  <f,\p>  (or  <(|),g>  ).  Thus,  by 
addition,  the  real  functional 

J('t',4>)  = R(V,4))  + R('f,4') 

* -<t,A*>  - <*,A*f>  ♦ <V,f>  + <f,'|'>  ♦ <g,(|»  ♦ <*,g>  (2.15) 

is  a bivariational  approximation  to 

o(«l^.4')  = <g.<)»  + <f,^>  = z(f,g;i;)  (2.I6) 

with  the  property 

J(y,«l»)  = Z(f,g;t)  - <6il;,A(5(|»  - <6^,A*6i|*>  . (2.17) 
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3.  MIXED  VARIATIONAL  BOUNDS  ON  Z 


In  the  event  that  H (the  self-adjoint  part  of  A)  is  a positive 
operator,  with  an  inverse  H it  is  possible  to  construct  two  special 
cases  of  the  bi variational  functional  J(4',*)  which  provide  complementary 
(upper  and  lower)  variational  bounds  on  Z(f,q;?)  . To  do  this,  we  think 
in  terms  of  the  'mixed'  vectors 

X = X(|)  + , 

y = X(J)  - pij;  , 

where  the  scalar  multipliers  X and  y are  real  and  satisfy 

2Xu  = 1 . (3.3) 


(3.1) 

(3.2) 


Combining  equations  (2.4)  and  (2.6),  we  see  that  x and  y are  the 
solution-vectors  of  the  simultaneous  equations 

(3.4) 

(3.5) 


Hx  + iu)y  = Xf  + wg  = p (say)  , 


and 


Hy  + iojx  = Xf  - ug  * q (say)  . 


If  we  now  choose  trial  vectors 
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X 


X4>  + u'i'  , 


(3.6) 


Y = A4 

- liH-  , 

(3.7) 

6x  = 

X - X = 

X64)  + 

(3.8) 

6y  = 

Y - y = 

A6(|)  - 

(3.9) 

we  find  that,  using  (3.3), 


J(H',<t) 


and 

I(X,Y) 


I(X.Y)  = -<X,HX>  + <Y,HY>  - <X,iioY>  + <Y,iuX> 

+ <X,p>  + <p,X>  - <Y,q>  - <q,Y> 

Z(f,g;C)  - <6x,(H6x+  iw6y)>  + <6y,(H6y+  iu6x)>  , 


(3.10) 


(3.11) 


Accordingly,  if  the  trial  vectors  X and  Y are  constrained  to 
satisfy  the  equation 


HX  + 1u)Y  * p (3.12) 

in  line  with  (3.4),  so  that 

H6X  + io)6y  = 0 , (3.13) 

it  follows  that  the  resulting  functional  I^(X,Y)  has  the  property 

I^(X,Y)  = Z(f,g;0  + <6y . (H+uV ^6y>  , (3.14) 

and  is  thus  a variational  upper  bound  on  Z . Similarly,  if 
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HY  + iojX 


(3.15) 


= q . 

in  line  with  (3.5),  so  that 

H6y  + io)6x  = 0 , 

it  follows  that  the  resulting  function  I_(X,Y)  has  the  property 

I (X,Y)  = Z(f,g;C)  - <6x, (H+oj2h" * )6x>  , 

and  is  thus  a variational  lower  bound  on  Z . 

In  terms  of  an  arbitrary  trial  vector  Y , we  have 

I^(X.Y)  = I(H'‘(p-iajY),Y)  = K^(Y)  say. 

wi  th 

K^.(Y)  = <Y,HY>  + <(p-ia)Y),H"^(p-iujY)>  - <Y,q>  - <q,Y>  . 

Similarly, 

I.(X,Y)  = I(X,H"^(q-iuX))  = K_(X)  say, 

with 

K_(x)  = -<X,HX>  - <(q-ia)X),H'^(q-iu)X)>  + <X,p>  + <p,X>  . 

Thus  we  obtain  the  complementary  upper  and  lower  variational  bounds 
K_(X)  < K_(x)  » Z(f,g;i:)  - K^(y)  < K^(Y)  . 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 
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with  the  bounding  functionals  each  depending  on  a single  'mixed'  complex 
vector.  Hence  we  call  them  'mixed'  variational  bounds.  It  is  a 
straightforward  matter  to  optimize  the  functionals  with  respect  to 
parameters  multiplying  the  trial  vectors,  and  if  H is  real  there  is 
separation  of  contributions  from  the  real  and  imaginary  parts  of  the  trial 
vectors.  The  ratio  X:p  is  also  a disposable  parameter. 

These  mixed  variational  bounds  hold  good  whenever 

>0.  n?«0,  (3.23) 

Eq  being  the  lowest  energy  eigenvalue  of  h . In  the  case  n = 0 , they 
hold  whenever 

El  - Eq  + V > 0 , (3.24) 

El  being  the  closest  eigenvalue  to  Eg  , since  Hq  does  not  contain  6g. 
If  Eg  is  degenerate,  then  Ei  = Eg  . Given  that  (3.23)  or  (3.24)  holds, 
we  have 


<4>,(h-E^+v)4»  = «I>,H(I»  > b<4-,4»  , b > 0 , (3.25) 

for  all  ♦ E , with 
n 

b = Eg  - E^  + V,  n ^ 0;  b = E,  - Eg  + v,  n = 0 . (3.26) 
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It  follows  from  (3.25)  that  H is  positive,  and  since  it  is  boundeo 
below  away  from  zero  the  inverse  operator  H ' exists  with  domain 
the  whole  of  . 

The  idea  of  introducing  a mixture  of  equations  like  (2.4)  and  (2.6) 
in  order  to  obtain  bounds  has  been  exploited  for  dissipative  systems  in 
real  spaces  by  Collins  [IJ.];  see  also  Herrera  [i4,lj^].  The  bi variational 
functional  I ( X.Y)  shows  the  saddle-type  dependence  on  X and  Y which  is 
necessary  for  complementary  bounds  [16,17]. 


4.  AVOIDANCE  OF  H~\  IMPLICIT  BIVARIATIONAL  BOUNDS 

A practical  disadvantage  of  the  bounding  functionals  K^(Y)  and 
K (X)  is  that  they  invol ve  the  inverse  , which  only  in  elementary 

cases  is  likely  to  have  a representation  simple  enough  to  permit  the 
evaluation  of  the  relevant  inner  products.  One  way  of  avoiding  H ^ in 
K^(Y)  is  to  write  as  a functional  of  X via  (3.12)  giving 


K^(Y(X))  = I(X,^HX-p))  = <(HX-p),H(HX-p)>  + <X,HX> 

- i <(HX-p),q>  + i <q,(HX-p)>  . 


(4.1) 


Simi larly. 


K.(X(Y))  - I(^HY-q),Y) 


- <(HY-q),H(HY-q)>  - <Y,HY> 


♦ i <(HY-q),p>  - ^ <p,(HY-q)>  . 


(4.2) 
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However,  these  forms  involve  (as  well  as  factors  of  w*'),  and  we  rule 
them  out  as  impractical. 

A better  way  of  avoiding  the  terms  in  (3.19)  and  (3.21)  is  to 

bound  them  separately,  using  individual  variational  bounds  of  the  type 

< -<X,Hx>  + <X.<^>  + «^.X>  + ;^l|Hx-i|F  . «-.XcH^,  b>0.  (4.3) 

which  follow  from  the  positivity  hypothesis  (3.25).  Taking  £ = p-iwY 
and  X = X in  (4.3),  we  find  that  (3.19)  gives,  after  simplification, 

K^(Y)  < I(X,Y)  + i||HX+iu)Y-p||2  . (4.4) 

Similarly,  putting  £ = q-itoX  and  x = Y (4.3),  we  obtain  from  (3.21) 
the  result 

K_(X)  > I(X,Y)  - il|HY+iu.X-q|P  . (4.6) 

Dropping  the  tildes  in  (4.4)  and  (4.5),  we  see  from  (3.22)  that  for 
arbitrary  vectors  X and  Y in  , 

I(X,Y)  - i||HY+iu)X-q||2  < Z(f,g;c)  < I (X,Y)  + il|HX+ia)Y-p||2  . (4.6) 

These  are  bivariational  bounds  on  Z(f,g;c)t  which  with  prescience  we  might 
have  derived  directly  from  (3.11)  and  (3.25).  They  hold  whenever  H is 
positive  and  bounded  below  away  from  zero,  so  that  a suitable  positive 
b can  be  found  according  to  (3.26). 
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The  mixed  vectors  (x,y)  and  (X,Y)  were  introduced  with  the  object  of 
deriving  variational  bounds  depending  on  a single  trial  vector,  and  so 
their  usefulness  has  evaporated  in  (4.6).  Referring  back  to  the  original 
vectors  and  (f.'t)  , the  bi variational  bounds  (4.6)  become 

J('f.-t)  - 1 {x2||A$-fl|2+y2||A*4'-g|P  - Re<A4>-f,A*'f-g>}  < Z(f,g;;) 

(4.7) 

1 J('f.<»)  + ^ {X2|lA4-fl|2+u2  1lA*T-g||2  + Re<A*-f .A*4'-g>}  . 

With  the  optimal . choice  for  the  ratio  X:p  of 

A|lA4>-f||  = p|lA*'f-g||  , (2Xy  = 1)  . (4.8) 


we  obtain  the  result 


J('f.«)  + ^ S(H',*)  - ^ C('f,4)  < Z(f,g;c)  < J(4',4)  + ^ 5(4', 4)  + ^ C(f,4)  , 


(4.9) 

where,  in  terms  of  H , 


J(4',«)  = -<y,H4»  - <(J»,H4'>  ♦ iai{<*,y>  - <y,4») 

(4.10) 

+ <y,f>  + <f,y>  + <g,^>  ♦ <4>,g>  , 

S(y,y)  ■ Real  <(»tHto)*-f,(H-1u>)y-g>  = Real  <A64»,A*6i)»>  , (4.11) 

I 

c(y,^)  = ||(H+Me-fl||(H-ia))y-g||  - 1M|  ||a*6.)»||  . (4.12) 
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We  call  the  bivariational  bounds  (4.9)  impi icit.  because  they  are 
contained  in  the  mixed  bounds  (3.22). 

In  the  special  case  of  zero  w , when  A becomes  self-adjoint,  the 
implicit  bounds  (4.9)  are  actually  tighter  than  others  previously  derived 
for  self-adjoint  operators  in  real  spaces  [l^,  see  also  19]- 

5.  EXPLICIT  BIVARIATIONAL  BOUNDS 

The  bounds  in  (4.9)  only  hold  when  the  operator  H is  positive  and 
bounded  below  away  from  zero  by  a positive  b , given  by  (3.26).  However, 
irrespectively  of  whether  this  condition  is  met,  it  is  possible  to  derive 
explicitly  alternative  bivariational  bounds  which  merely  require  the 
condition 


II A* II  > a||<l)||  , a > 0 , for  all  4>  e . 


Since 


(5.1) 


llA't'lp  = <A4>,A*>  = «t,A*A4»  = + u)2«I),4» 

1 (Ep.-V^)"  II^IF  I^IP. 


(5.2) 


where  E^,  is  the  energy  eigenvalue  of  h which  minimizes  (E^,-E^+v)^  , 
we  see  that  condition  (5.1)  is  satisfied  by  taking 


a^  * , a > 0 . 


(5.3) 
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We  bear  in  mind  that  E , ^ (unless  E„  is  deqenerate),  since  H does 

n n n n 

not  contain  6^  . Thus,  disregarding  the  exceptional  case  of  zero  ; and 
degenerate  E^  , we  can  always  find  a constant  a to  satisfy  (5. 1 ) . 

Applying  (5.1)  to  (4.12),  and  then  using  Schwarz's  inequality,  we  have 

C(f,4>)  > alMl  ||A6i|»||  > a|<64»,A*6ij;>l  . (6.4) 

The  magnitude  of  the  complex  number  on  the  right  of  (5.4)  is  greater  than 
or  equal  to  the  magnitude  of  its  real  part,  and  so  from  (2.17)  it  follows 
that 


C(4',4>)  > y2a|<6(|),A*6i|j>  + <6i|;,A6(j)>|  = V2a IZ-J('»',4>)  | . (5.5) 


Rearranging  (5.5)  we  obtain  at  once  the  explicit  bivariational  bounds 


J('f,<l>)  - 7 C(f,*)  < Z(f,g;^)  < 0(4',4>)  + J C(y.O) 

d ““  d 


(5.6) 


Existence  theorems  for  bounds  of  this  type  have  recently  been  presented, 
together  with  some  applications  [20,21]. 

6.  THE  CASE  f = g 

When  f = g,  as  is  the  case  for  the  dynamic  polarizabil it> 
in  (1.1),  the  equations  (2.4)  and  (2.6)  become 


(H+i(j)0  » f = (H-iu))i^  . 


(6.1) 
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If  d similar  relationship  is  imposed  on  the  trial  vectors  and  'll  by 
taking 


♦ = (H-iu)0,  f = (H+ioj)0  , (6.2) 

for  some  trial  vector  0 which  is  supposed  to  approximate  *f  , 

the  functionals  in  (4.9)  and  (5.6)  become 

J = -2<0,H(H2+uj2)0>  + 2<f,HO>  + 2<H0,f>  , (6.3) 

and 

C = S = |l(H2+a)2)0.f||2  . (6.4) 

The  bi variational  bounds  (4.9)  and  (5.6)  become  variational  bounds, 
involving  high  powers  of  H , of  a type  previously  obtained  [3,5].  The 
and  H**  rule  them  out  for  practical  purposes.  Thus  there  seems  little 
point  in  trying  to  impose  a constraint  like  (6.2).  However,  when  choosing 
trial  vectors,  it  would  be  sensible  to  have  in  mind  the  relationship 
between  the  respective  w-dependence  of  ♦ and  'f  which  is  implied  by 
(6.2). 

There  are  some  simplifications  to  be  made  if  f is  a real  vector, 
and  the  Hamiltonian  h (and  hence  H ) is  a real  operator.  Then  it 
follows  from  (6.1)  that  » 0 . Accordingly,  let  us  set 

♦ * + i*2»  y = - 1*2  * ♦ » (6.5) 
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in  the  functionals  J,  C and  S .where  and  4>2  are  real  vectors. 


We  obtain  the 

functionals 

J(4',4>) 

= -2<$i,H^i>  + 2<^2*H^2^  4u)<4>j,$2^  + 4<4>i,f>  , 

(6.6) 

C($,4>) 

= |lH4>i-u)4>2-f  ||2  + II H$2'*’<*>4>i  IP  , 

(6.7) 

and 

S(4>,4>) 

= ||H4>i-a>4>2-f  Ip  - ||H'I>2+<jD'1>i  IP  . 

(6.8) 

The  implicit  bivariational  bounds  (4.9)  become 

J(4>,4>)  - ||^i4>2+(D4>i|p  < Z(f,f;c)  < + ||^^4)l-(iJ4>2-f||2  , (6.9) 

and  the  explicit  bounds  (5.6)  take  the  form 

j ($,$).  i C(*,4>)  < Z(f,f;c)  < J(4>,<1>)  + I C(5,<t')  (6.10) 

which  at  v * 0 is  essentially  that  given  by  Burrows  [8].  The  mixed 

bounds  K^(Y)  and  K_(X)  simplify  in  this  situation,  too,  particularly 

_ V 

if  we  take  X = y = 2 ^2  ^ so  that  X = , Y = i4>2^  , p = fv7 


and  q = 0 , yielding 

K^(Y)  « 

2<4>2,H*2>  + 2 <(f+u)4>2).H"*(f+w4>2)> 

(6.11) 

and 

K.(X)  * 

-2<*j,H^i>  - 2u)2<4>j,H  *4>j>  + 4<4>j,f>  . 

(6.12) 
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Goscinski  [^]  has  qiven  the  amplitude-optimized  version  of  the  bound  in 
(6.12)  for  V = 0 . 


7.  DISCUSSION 

Interest  in  the  mixed  variational  bounds  K^(Y)  and  K_(X)  is 
primarily  theoretical.  Not  only  do  they  contain  the  implicit  bivariational 
bounds,  but  also  they  can  lead  to  bounds  in  terms  of  other  known  quantities 
To  take  a simple  example,  when  ui  is  small  the  solution  to  equation  (2.4) 
is  approximately  H"*f  - iuH'^  f . Thus  if  we  take 


♦ i = CiH’V, 


♦z 


-iojc 


f 


(7.1) 


in  the  simplified  functionals  (6.11)  and  (6.12),  and  optimize  with  respect 
to  Cl  and  C2,  we  should  obtain  bounds  on  Z(f,f;c)  which  are  accurate  for 
small  w . The  bounds  are 


(7.2) 


in  terms  of  the  'sum  rules' 


Sn  = <f,H"'^V>  . 


(7.3) 


For  ground-state  hydrogen,  the  bounds  (7.2)  give  the  result  (in  atomic  units) 


4.2490  < a(iui)  < 4.2503 


(7.4) 


for  the  dipole  polarizabi 1 ity  at  v » 0,  w * 0.1. 
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When  H 1s  positive  the  quantity  Z(f,f;c)  Is  a serles-of-Stleltjes- 
representable  function  of  , and  with  suitable  choices  of  trial  vector 
different  families  of  Pad^  approxlmant  bounds  can  be  derived  from 
and  K_  [6,^].  More  generally,  whenever  <f,(H2+u2)’*g>  Is  real, 
it  can  be  shown  that 

Z(f.g;c)  = <(H%).(H2+tu2)"MH%)>  - <(H%).(H2+u)2)‘'(H'^^q)>  . (7.5) 

whence  Z(f,g;c)is  the  difference  of  two  series-of-Stiel tjes-representable 
functions.  The  and  K_  functionals  are  again  the  appropriate  ones  to 
yield  the  Pad6  approxlmant  bounds  for  this  kind  of  situation  [^],  rather 
than  bivariatlonal  functionals  which  lead  to  Pad^  approximants  plus 
correction  terms  [21 .24]. 

From  a practical  standpoint,  it  is  unlikely  that  a convenient 
representation  of  will  be  available,  in  which  case  the  mixed  bounds 
K_j_(Y)  and  K_(X)  are  not  of  direct  interest  even  though  they  only  involve 
a single  trial  vector.  Likewise  bounding  functionals  containing  and 
higher  powers  can  be  ignored,  because  of  the  consequent  difficulties  in 
evaluating  the  inner  products  when  sensible  trial  vectors  are  employed. 

Thus  for  a practical  tool  we  are  left  with  the  bivariatlonal  bounds,  in 
either  the  implicit  form  (4.9)  (only  valid  when  H is  bounded  below  away 
from  zero  by  a positive  number  b)  or  the  explicit  form  (5.6)  (always 
valid).  Although  an  extra  trial  vector  is  involved,  only  H and  appear 
in  the  inner  products  and  this  advantage  is  crucial.  The  bivariatlonal 
bounds  should  still  be  used  even  when  f = g;  the  simpler  versions  (6.9) 
or  (6.10)  are  relevant  when  f * g is  a real  vector  and  H is  real. 
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Even  when  the  implicit  bounds  (4.9)  are  available,  they  will  not 
necessarily  yield  better  results  than  the  explicit  bounds  (5.6).  For 
example,  when  cj  is  large,  the  number  a (given  by  (5.3))  is  of  order 
0)  , whereas  the  number  b (given  by  (3.26))  is  not.  With  the  correct 
asymptotic  choices 

^ = zli  , ^ = il  (7.6) 

0>  0) 


the  explicit  bounds  give 

J tjc  = ^<f,g>  - <g,f>)  + ,-77(<f.Hq>  + <q,Hf>)  ± 0(,-^)  , (7.7) 

o OJ  uJ  tu 

_ 2 

whereas  the  bounds  (4.9)  leave  some  uncertainty  in  the  w term.  However, 
in  cases  like  (6.4)  when  (S-C)  is  zero,  or  very  small,  it  is  clear  that  the 
lower  bound  in  (4.9)  is  better  than  the  lower  bound  in  (5.6).  We  notice 
also  that  the  explicit  bounds  can  never  have  quite  the  correct  w-dependence 
because  of  the  square-root  defining  the  number  a . Thus,  when  both  are 
available,  the  implicit  and  the  explicit  bivariational  bounds  should  each 
be  investigated  in  any  given  situation  to  see  which  gives  the  better 
resul ts. 

Extensive  calculations  of  bounds  on  dynamic  polarizabilities  (f=g)  for 
two-electron  atoms  at  zero  w or  zero  v have  recently  been  carried  out  by 
Glover  and  Weinhold  [^-^].  The  functionals  of  Braun  and  Rebane  [2]  which 
they  employed  led  to  inaccuracies  for  small  values  of  w.  Applications  of  the 
Glover-Weinhold  techniques  to  the  bivariational  bounds  developed  here  are  being 
explored.  A pleasing  feature  of  the  present  approach  is  that  the  same 
bivariational  functionals  provide  bounds  whether  or  not  tij=0,  v=0,  or  f=g. 
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